Mass of the Be Meson in Three-Flavor Lattice QCD 



Ian F. Allison,^ Christine T.H. Davies,^ Alan Gray,^ Andreas S. Kronfeld,'^ Paul B. Mackenzie,"^ and James N. Simone'^ 

(HPQCD, Fermilab Lattice, and UKQCD Collaborations) 

^Department of Physics and Astronomy, Glasgow University, Glasgow, Scotland, United Kingdom 

^Department of Physics, The Ohio State University, Golumbus, Ohio, USA 
''Theoretical Physics Department, Fermi National Accelerator Laboratory, Batavia, Illinois, USA 

(Dated: November 19, 2004) 

We use lattice QCD to predict the mass of the Be meson. We use the MILC Collaboration's 
ensembles of lattice gauge fields, which have a quark sea with two flavors much lighter than a third. 
Our final result is ms^ = 6304 ± 12t^^ MeV. The first error bar is a sum in quadrature of statistical 
and systematic uncertainties, and the second is an estimate of heavy-quark discretization effects. 
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Recently there has been a significant breakthrough in 
numerical lattice calculations of QCD With new, 
improved techniques for incorporating light sea quarks, 
lattice QCD agrees with experiment at the few percent 
level for a wide variety of quantities. This progress sug- 
gests that lattice QCD could play a big role in particle 
physics, especially as an aid to understanding the flavor 
sector of the Standard Model . 

In flavor physics, the central aim is to search for evi- 
dence of new phenomena. Before applying results from 
numerical lattice QCD for such purposes, it is helpful to 
have as many tests as possible. Although lattice gauge 
theory has a solid mathematical foundation, numerical 
simulations are not simple. The impressive results of 
Ref. |lj have been achieved only with the fastest method 
for simulating light quarks. The price for speed is an 
unproven assumption (discussed below), which clearly 
warrants further scrutiny. In addition, the cutoff effects 
of heavy quarks are controlled using effective field theo- 
ries. Although most heavy-quark phenomenology relies 
on this framework, it is important to find out how well 
it describes discretization errors in lattice calculations. 

The ideal way to test a theoretical technique is to 
predict a mass or decay rate that is not well-measured 
experimentally, but will be measured precisely soon. 
Some examples are in leptonic and semileptonic decays 
of charmed mesons, which are being measured in the 
CLEO-c experiment. They are sensitive to both the light- 
quark and heavy-quark methods, and are under investi- 
gation 

Another example, pursued here, is the mass of the 
pseudoscalar B^. meson, the lowest-lying bound state of 
a bottom anti-quark (6) and a charmed quark (c). The 
Be mass principally tests the heavy-quark methods of 
lattice QCD. Based on experience with hh 5J and cc 
mass splittings, we expect only mild sensitivity to the 
light quark mass (of the sea quarks) once the mass is 
small enough to allow uninhibited creation and annihila- 
tion of virtual light quark pairs. Preliminary versions of 
this work have been given at conferences 



Until now. Be has been observed only in the semilep- 
tonic decay B^ — > Jjjjjl~^vi, with a mass resolution of 
around 400 MeV d During Run 2 of the Fermilab 
Tevatron, Be is expected be observed in non-leptonic de- 
cays, with a mass resolution estimated to be 20-50 MeV 
llOl . Our total uncertainty is much smaller than the cur- 
rent experimental accuracy, and comparable to the pro- 
jections, so we may claim to be predicting the mass of 
the Be meson. 

Heavy-quark discretization effects are a challenge, be- 
cause feasible lattice spacings a are about the same as 
the Compton wavelength of the bottom and charmed 
quarks. The distances are both shorter than the typical 
distance of QCD, which is about 1 fm. The obvious strat- 
egy is to use effective field theories to separate long- and 
short-distance scales. This reasoning has led to the de- 
velopment of non-relativistic QCD (NRQCD) for quarko- 
nium and heavy-quark effective theory (HQET) for 
heavy- light mesons |l2l | . In lattice gauge theory, this rea- 
soning has led to two systematic methods for discretizing 
the heavy-quark Lagrangian: lattice NRQCD Jd,,J^] and 
the Fermilab heavy-quark method Il5| . A strength of 
both is that the free parameters of the lattice Lagrangian 
can be fixed with quarkonium. Then, with no free param- 
eters, one obtains results for heavy-light systems (such as 
D and B mesons). The same procedure applies here: we 
obtain m with the same bare quark masses that repro- 
duce the bottomonium |^ and charmonium spectra. 

It is beyond the scope of this Letter to review the de- 
tails of heavy quarks in lattice gauge theory .16J. T he 
couphngs of the Lagrangian are adjusted so that [l^ 

Aat = Cqcb + Smiji^h^ + h^h^) + 

Y,a'"UmQa)On (1) 

n 

where = can be read "has the same mass spectrum as." 
The 5m term is an unimportant overall shift in the mass 
spectrum; /i+ {h~) is a effective field for quarks (anti- 
quarks); the On are the effective operators of the heavy- 
quark expansion, of dimension dimO„ = 4 -f s„, s„ > 1; 
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and a is the lattice spacing. The coefficients /„ arise from 
the short-distance mismatch between lattice gauge the- 
ory and continuum QCD. By choosing an improved lat- 
tice Lagrangian £iat, the can be reduced. In practice, 
however, one must vary a and also estimate the effects of 
the leading 0„ on the mass spectrum. 

Our calculation employs an idea from a quenched cal- 
culation ^3 (omitting sea quarks), namely to use lattice 
NRQCD for the b quark and the Fermilab method for 
the c quark. The lattice NRQCD Lagrangian has a 
better treatment of interactions of order u^, where v is 
the heavy-quark velocity. The Fermilab Lagrangian 0| 
has a better treatment of higher relativistic corrections, 
which is helpful since the velocity of the c quark in Be 
is not especially small, v'^ « 0.5. Thus, we expect this 
combination to control discretization effects well. This 
choice also means that our calculation directly tests the 
heavy-quark Lagrangians used in Ref. 0. 

We work with ensembles of lattice gauge fields from the 
MILC Collaboration . Each ensemble contains several 
hundred lattice gauge fields, so statistical errors are a few 
per cent. The gluon fields interact with a sea of "2 -I- 1" 
quarks: one with mass tuned close to that of the 
strange quark, and the other two as light as possible. In 
this work we use ensembles with light mass mi = O.Itos, 
mi = 0.2ms, and m; = 0.4ms. The gluon and sea-quark 
Lagrangians are improved to reduce discretization effects. 
We use three lattice spacings, a ^ jj, ^, fm. Further 
details are in the MILC Collaboration's papers [l^ . 

A drawback of the MILC ensembles is that the sea 
quarks are incorporated with "staggered" quarks. A sin- 
gle staggered quark field leads to four species, or "tastes," 
in the continuum limit. Sea quarks are represented (as 
usual) by the determinant of the staggered discretization 
of the Dirac operator. To simulate 2 tastes (1 taste), the 
square root (fourth root) of the 4-taste determinant is 
taken. The validity of this procedure is not yet proven 
for lattice QCD, although a proof does go through in at 
least one (non-trivial) context Moreover, one finds 
that interacting improved staggered fields split into quar- 
tets lio], as is necessary. Since our prediction of the Be 
mass tests this ingredient of the calculation (albeit indi- 
rectly), we do not assign a numerical error bar to this 
issue. 

As in Ref. 0|, we calculate mass splittings, namely 

A^T = ms, - {m^ + mr)/2, (2) 

Ad,b, = ms, ~ {mo, + ffiBj, (3) 

where m^ = (m,,^ -|-3mj/^)/4, mn, = {mo, -l-3mD*)/4, 
and ifiBs = (jfiBs + 3mB*)/4 are spin-averaged masses. 
We refer to (m^ -I- mx)/2 and {mjj^ + friB^) as the 
"quarkonium" and "heavy-light" baselines, respectively. 
Our result for ms^ comes from our calculated oA^x and 
aAjj^Bs (in lattice units), combined with the lattice spac- 
ing a and the experimental measurements of the base- 
lines. We use the 2S-1S splitting of bottomonium to 



define a, but on the MILC ensembles several other ob- 
servables would serve equally well Fll. 

Many uncertainties cancel in mass splittings. Lattice 
calculations integrate the QCD functional integral with 
a Monte Carlo method, and the ensuing statistical er- 
ror largely cancels when forming a difference. The mass 
shifts Sm in Eq. drop out. The spin-averaging cancels 
the contribution of the hyperfine operator h^iH ■ Bh^. 
(We do not spin-average T with rji,, because the latter 
remains unobserved.) The discretization errors from fur- 
ther terms in Eq. cancel to some extent, especially 
with the quarkonium baseline. Most crucially, all masses 
in Eqs. Q and (O are "gold-plated" 0, in the sense 
that the hadrons are stable and not especially sensitive 
to light quarks. (Hence we use Dg and Bg, not D and B.) 

We turn now to a discussion of our numerical work. 
First we discuss briefiy how to compute the meson 
masses. Then we consider systematic effects that can 
be addressed directly by varying the bare quark masses 
(light and heavy) . Finally, we consider the remaining dis- 
cretization effects, by changing the lattice spacing and by 
studying the corrections in Eq. Q). 

In lattice QCD, each meson mass is extracted from a 
two-point correlation function, which contains contribu- 
tions from the desired state and its radial excitations. 
We use constrained curve fitting [ij, usually including 
5 states, but checking the results with 2-8 states in the 
fit. We find that the extraction of the raw masses is 
straightforward on every ensemble. 

Statistical errors are obtained with the bootstrap 
method. The statistical precision on A^x is about 4% 
and on Ad,b, about 1.5%. But since A^x ~ 40 MeV 
and Ad^Bs ~ —1200 MeV, the statistical error on m^^ 
ends up being much larger with the heavy-light baseline. 

Figure n shows how the splittings depend on the light 
quark mass m;, for the ensembles with a ~ ^ fm. The 
dependence on m; is hardly significant. We extrapolate 
linearly in mi/mg, down to the value that reproduces the 
pion mass 0|- The mild dependence on m; also suggests 
that the uncertainty from the known (but small) mis- 
tuning of the strange quark sea is completely negligible. 

The bare masses of the heavy quarks are chosen as fol- 
lows. Since the overall mass is shifted [by Sm in Eq. (^)], 
we compute the kinetic energy of bb and cc mesons of 
(small) momentum p, and choose the bare b and c quark 
masses so that it is /2m, where m is the physical QQ 
mass. The statistical and systematic uncertainties of the 
kinetic energy imply a range of bare quark masses. We 
compute the effect on B^ for different bare b and c masses 
and derive an error of 10 MeV (5 MeV) in A^x and 
^DsBs from this source. 

Figure 121 shows how A^x depends on lattice spacing a. 
The change is insignificant. Lattice spacing dependence 
stems from all parts of the lattice QCD Lagrangian. In 
our case, the heavy-quark discretization effects, espe- 
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FIG. 1: Sea-quark mass dependence of A^x and AdsBs- 



cially for the c quark, are expected to dominate. Un- 
fortunately, the dependence on rrica [of the coefficients in 
Eq. Q] does not provide a simple Ansatz for extrapola- 
tion. 

We shall treat discretization errors with Eq. Q, us- 
ing calculations of the short-distance mismatch and esti- 
mates of the On- This approach is itself uncertain, but 
it is preferable to ignoring the issue. The results of such 
an analysis are in given in Table and the following 
paragraphs explain how the entries are obtained. 

As usual, we classify the operators On in Eq. ^ ac- 
cording to the power-counting scheme of NRQCD (or, for 
Ds and mesons, HQET). Table n lists those of order 
in NRQCD; in HQET they are of order l/m^, n = 
1,2,3,3. The spin-orbit interaction /i^^zS • {D x 
is omitted, because its matrix elements vanish in the S- 
wave states considered here. 

The contribution of the hyperfine interaction /i^iS ■ 
Bh^ cancels for spin-averaged masses to, by construc- 
tion, but we must still estimate its effect on tox and TOb^. 
In the heavy-quark Lagrangians we are using, the hy- 
perfine coupling is correctly adjusted only at the tree 
level. Indeed we find discrepancies in the hyperfine split- 
tings TOfl. — m]j_. and rrij/^ — for the c quark and 
tob* — niB^ for the b quark. The size of the discrepancy 
agrees with the expectation from the one-loop mismatch 

r , 1 , , ^ , 1 1 , ^ , , 1 1 ^ , , 1 1 q 
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FIG. 2: Lattice-spacing dependence of A^t. 



in the coefficient. The hyperfine entries for tox and ms^ 
are obtained by combining the coefficient mismatch with 
the computed hyperfine splittings. 

For TOSej and ^tox, the matrix elements of the 

Darwin term h^D ■ Eh^ and the relativistic corrections 
(-D^) and h^Dfh^ are obtained from po- 

tential models. For ffiOs and fhs^ we use HQET di- 
mensional analysis: {D ■ E) ^ A"^, (Z3^) ^ A^, with 
A = 700 MeV. Next we multiply the estimated matrix 
elements (On) with the mismatch coefficients fnimga). 
We have explicit tree-level calculations of them for the 
Fermilab Lagrangian used for the c quark. For the b 
quark the mismatch starts at order a^, so we take /„ to 
be of order as with unknown sign. The resulting shifts 
from the c quark are larger, but their sign is definite. 

The entries in Table U for (1)2)2 jji 

are uncer- 
tain. The cancellations across each row are reliable, 
but the overall magnitude could be larger. The same 
potential model suggests a shift in our mh^ — of 
about —10 MeV, consistent with the computed discrep- 
ancy Thus, the charmonium spectrum suggests 
that the entries are reasonable. 

Table m suggests that our results for to_b^ will be too 
low, and that ttib^ will be lower with the heavy-light 
baseline than with the quarkonium baseline. We could 
apply the shifts in Table to our lattice QCD results. 
Our aim, however, is to test lattice QCD. Therefore, we 
treat these shifts not as corrections but as uncertainties. 
Since we claim to know the sign in the important cases, 
the associated error bars are asymmetric. Repeating this 
analysis at other lattice spacings yields consistent error 
estimates. 

After extrapolating the light quark mass and accu- 
mulating the other systematic uncertainties we find (at 
a — fm) 

A^x = 39.8±3.8± 11.21^^ MeV, (4) 
Ad^b, = - [1238±30± Uig^] MeV, (5) 

where the uncertainties are, respectively, from statistics 
(after extrapolating in mi/ms), tuning of the heavy- 
quark masses, and heavy-quark discretization effects. 
The results for A^x at a = jj, ^ fm are completely 



TABLE I: Estimated shifts (in MeV) of masses and splittings 
A^T and AogBa at a = ^ fm. Entries show what should 
be added to the masses and splittings to compensate for dis- 
cretization errors. Dots (■ ■ ■) imply the entry is negligible. 
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FIG. 3: Comparison of theoretical work, with references in 
brackets and our equation numbers in parentheses. 



consistent. For the Be mass we find 

ms, = 6304 ± 4±llli^MeV, (6) 
me, = 6243±30±lll^^ MeV, (7) 

restoring, respectively, the quarkonium and heavy-quark 
baselines. We have carried out more checks on the quark- 
onium baseline, so we take Eq. @ as our main result. 
Given the rough nature of the last error bar, we consider 
the agreement of the two results to be reasonable. Fur- 
ther work with more highly improved Lagrangians and 
at finer lattice spacing should reduce this error. 

Our results are compared to other theoretical pre- 
dictions in Fig. 121 including potential models Esl . 
quenched lattice QCD and potential NRQCD 
I25I |26| |. The quarkonium baseline is shown for refer- 
ence. Our result is so much more accurate than the pre- 
vious lattice QCD result simply because we have 
eliminated the quenched approximation. If our predic- 
tion, Eqs. lO and |(7J), is borne out by measurements, 
it lends confidence in lattice QCD, not only in MILC's 
method for including sea quarks, but also in the control 
of heavy-quark discretization effects using effective field 
theory ideas. Moreover, within this framework it is clear 
how to improve the lattice QCD Lagrangian to reduce 
the remaining uncertainties. 
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